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Abstract. We present a model for the interpretation of 
the 7-ray background emission from the Galactic halo, 
which we attribute to interaction of high-energy cosmic 
rays with dense molecular clumps. In a wide range of 
clump parameters we calculate the expected 7-photon 
flux, taking into account for the first time optical depth 
effects. This allows us to derive new constraints on masses 
and sizes of possible molecular clumps and their contribu- 
tion to the galactic surface density. 

The observed 7-ray distribution can be explained best 
by models with a flattened halo distribution of axis ratio 
~ 0.3. If optical depth effects are important, the clumps 
must have radii of ~ 6 AU and masses of ~ 1O _3 M0. 
This would result in a total mass of ~ 2 • lO n M0 for such 
clumps and contribute with X ~ 140 M Q pc -2 to the local 
surface density. 
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1. Introduction 

Recent observational results of the anisotropy of the cos- 
mic microwave background, deuterium abundance from 
cosmological nucleosynthesis, dynamics of clusters of 
galaxies and the data from the Supernova Cosmology 
Project provide evidence for the distribution of mass in 
the universe in the proportion: f2b — 0.05 for baryons, 
tt m ~ 0.35 for non-baryonic matter, and £Ia — 0.6 for 
the cosmological constant (Turner, 1999). Since luminous 
baryons contain only f2£ — 0.007, more than 85 % of the 
baryonic mass is still undetected. Pfenniger & Combes 
(1994) and more recently Walker & Wardle (1998) have 
suggested that a considerable fraction of baryonic mass 
can be contained in dense molecular clumps of AU sizes 
and about a Jovian mass, which in turn may manifest 
themselves through extreme scattering events (ESEs) - 



dramatic flux changes of compact radio quasars over sev- 
eral weeks (Fiedler et al., 1987). However, estimates of 
this fraction are uncertain because the existing statis- 
tics of the ESEs sets a wide range for the covering fac- 
tor of the radio-refracting regions: / ~ 10~ 4 to 5 ■ 10 -3 
(Walker & Wardle, 1998). From this point of view 7- 
rays produced via interaction of high energy cosmic rays 
(CR) with nucleons of dense clumps through the process 
p+p^p + p + n° p + p + 2j are a unique probe of 
baryons hidden in dense H2 clumps. 

Recently several groups have studied the 7-ray emis- 
sion from H 2 clumps in detail using EGRET [E > 100 
MeV photons) data (Chary & Wright 1999, De Paolis et 
al. 1999, Dixon et al. 1998, and Sciama 1999). An essential 
assumption of all these attempts is that the dark matter in 
the halo is optically thin to both, the exposing high-energy 
protons, and the resulting 7-ray photons. However, as we 
will argue below, most of the possible baryonic dark mat- 
ter candidates are dense and compact enough to allow 7- 
ray emission only from thin external skin layers, and thus 
the existing EGRET data can trace a small fraction /t of 
baryons in the halo only. In this Letter we determine the 
7-ray background emission from small optically thick H2 
clouds. 

2. 7-ray emission from dense clumps 

At present, no direct observations of cold, T ~ 3K, dense 
molecular clumps which might carry a considerable frac- 
tion of baryonic dark matter do exist, and their parame- 
ters, such as masses and radii, may vary in a wide range. 
Pfenniger & Combes (1994) have argued that a consider- 
able fraction of baryonic dark matter in galaxies can be 
accounted for by dense clumps of predominantly molecu- 
lar hydrogen of Jovian masses and radii of 30 AU. Ger- 
hard & Silk (1996) estimated masses of ~ 1 M Q and radii 
of - 0.03 - 0.1 pc. Walker & Wardle (1998) have shown 
that the extreme scattering events might be naturally ex- 
plained if the clumps have radii of ~ 3 AU. Larger radii, 
~ 10 AU, are proposed by Draine (1998) from consid- 



2 



P.M.W. Kalberla, Yu.A. Shchekinov, and R.-J. Dettmar: H2 dark matter in the galactic halo from EGRET 



eration of optical lensing of stars by dense gas clouds. 
Quite different arguments, based on the analysis of turbu- 
lent motions of HI gas in the halo, lead Kalberla & Kerp 
(1998, KK98) to conclude that a significant fraction of the 
dark matter in the halo can be in form of gas clouds with 
masses of £ 2 x 10~ 3 M and radii of <~ 10 AU. To distin- 
guish between these different proposals, we need to check 
whether optical depth effects might affect the observed 7- 
ray emission. High energy (> 100 MeV) CR protons in 
H 2 are attenuated by a factor e for Scr — 40 g cm~ 2 
(Salati et al., 1996). The mass column density of a clump 
along the radius is iV c = 4.2 x 10 6 M/r 2 g cm~ 2 , where 
r c is the clump radius in AU, M, its mass in M Q . Using 
the parameters as proposed by Walker & Wardle (1998) 
we find from this crude estimate that N c exceeds Scr by 
two orders of magnitude. 7-photons produced inside the 
clumps suffer from absorption. According to Salati et al. 
(1996), an optical depth of one in H 2 is reached at S p ~ 80 
g cm -2 , twice as big as Scr- 

After this first estimate we define /t as the fraction of 
7-ray emission from a dense clump relative to the emission 
from an optical thin cloud. We derive /t by integrating 
the 7-ray emission from the sphere assuming that the in- 
falling cosmic rays are distributed isotropically. This in 
turn results in an isotropic distribution for the 7r°-photons. 
We assume that the density p of the clump is constant and 
derive 
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where Li = Si/ p, i = CR, p, 

_ e -(rc+r)/£«[ 1 + ( rc +,.)/£.]. 



(2) 



/t is an approximation only, since we use frequency av- 
eraged attenuation lengths Scr and S p , we neglect also 
electron Bremsstrahlung. 

In Fig. 1 we display /t as a function of the clump ra- 
dius r c for clumps with masses 10~ 4 to 10 _2 M©. There is 
significant absorption for radii r c < 30 AU. Clouds de- 
scribed by Pfenniger & Combes (1994) and by Gerhard & 
Silk (1996) are transparent and therefore in these models 
the 7-ray intensity is proportional to the total mass con- 
tained in dense clouds and clumps. On the contrary, dense 
clumps in models described by Walker & Wardle (1998), 
Drainc (1998), and KK98 are optically thick. In this case 
the determination of the mass of a baryonic dark matter 
halo from the observed 7-ray emission depends on /p. 

The 7-ray emission is caused by nuclear interactions 
between cosmic rays and matter (most prominent the n°- 
decay), by electron Bremsstrahlung, and by inverse Comp- 
ton interactions. We used the source functions for nucleon- 
nucleon interactions and Bremsstrahlung as published by 
Bcrtsch ct al. (1993). For the inverse Compton interaction 
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Fig. 1. /t, the fraction of the 7-ray emission from a clump 
with radius r c relative to the emission from an optical thin H2 
cloud. The curves are for clump masses 10~ 4 to 10~ 2 Mq. 



we have used the "galprop" database according to Strong 
& Moskalenko (1997). The total 7-ray flux observed on 
Earth is a superposition of fluxes produced by dark mat- 
ter halo clouds and diffuse components of the interstellar 
gas: the extra-planar diffuse ionized gas (Dettmar, 1992), 
the H 1 disk, and the gaseous halo with H 1 gas and plasma 
(KK98). The density Pi{R g ,z), (R g = ^Jx 2 +y 2 ), of the 
various gaseous components is given by 



Pi(R g ,z) = n i .gi( J R g )exp 



-$(*) 



af (1 + a l + P l ) 



(3) 



with rii, the local midplane density of the individual com- 
ponent, <fr(z) the gravitational potential, Oi the corre- 
sponding velocity dispersion, and and Pi according to 
Parker (1966) quotients which determine the pressure of 
the magnetic field and cosmic rays relative to the gas pres- 
sure. According to KK98 all gaseous components have a 
common radial density distribution according to 



9l (R g ) =sech 2 (R g /A 1 ) /sech 2 (R Q jA x ) , 



(4) 



with a radial scale length A\ = 15 kpc. This relation 
was introduced by Taylor & Cordes (1993) to describe 
the diffuse ionized gas component. The scale height of the 
gaseous halo according to KK98 is h z = 4.4 kpc. 

A widely used standard expression to describe the den- 
sity distribution in the galactic halo is (e.g. De Paolis et 
al., 1999) 



pn 2 (x,y,z) = p (q) 



R l 



a 2 + x 2 + y 2 + (z/q) 2 ' 



(5) 



where x, y, z are the galactocentric coordinates, i?o is the 
Solar distance, po(q), the local dark matter density. 
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3. Model calculations 

We discuss here several models for the distribution of dif- 
fuse gas components and the dark matter clumps in detail: 

1) the disk, the extra-planar diffuse ionized gas and the 
(observed) gaseous halo with parameters as in KK98 
model according to Eq. (|J) and (||). The local mid- 
plane density of the 7-ray emitting gas is rig = 0.0025 
cm -3 , and the CR density is derived from the pressure 
equilibrium with the observed gas. In this model there 
is no significant 7-ray emission from the halo. 

2) same as the previous model with the best fit local mid- 
plane density of the halo 7-ray gas rig = 0.065 cm~ 3 . 
In this model the flat rotation curve demands a local 
total midplane density rig — 0.7 cm~ 3 (KK98). As- 
suming that the total dark matter in the Milky Way 
is contained in H2 clumps, we obtain /x = 0.09. 
the dark matter halo with a distribution in the form 
(||) with a = 5.6 kpc, q = 0.3, R = 8.5 kpc, and 
a constant CR energy density in the halo of 0.12 eV 
cm -3 . This model is similar to the model of De Pao- 
lis et al., however without a central hole in the dark 
matter distribution at R < 10 kpc. We fit a local mid- 
plane density rig = 0.55 cm -3 . For a rotation velocity 



3) 



of Vq = 220 kms 



we derive rig 



0.5 cm 3 hence 



n o/ n o = l-l- According to Eq. (1) /x < 1, and within 
the errors we obtain for this model /x = 1. 

4) same as the previous model but with q = 1 and ng = 
0.18 cm -3 , we derive rig = 0.24 cm~ 3 corresponding 
to / T = 0.75. 

5) same as the model 4, but with a central hole for R < 
10 kpc. We determine rig = 0.35 cm -3 . This is the 
model proposed by De Paolis et al. (1999) with rig = 
0.32 cm~ 3 . This model cannot reproduce the rotation 
velocity observed in the inner galaxy; /x is undefined. 

We estimate the errors in the determination of nfj and rig 
to about 10%. The isotropic 7-ray background was fitted 



to 5(±1) • lO^ctss" 
the rate of 4(±1) • 10" 



cm sr , in good agreement with 
" 6 cts s -1 cm _2 sr -1 due to blazars as 



determined by Mukherjee & Chiang (1999). 



4. Discussion and conclusions 

For all models we calculated the RMS deviation between 
data and model, individually for each of the scans at con- 
stant longitude and further for all of the data under con- 
sideration (Fig. 3). We found that the RMS values derived 
this way were for all of the models biased by a number of 
positions with significant deviations between model and 
data. Predominantly the deviations were due to 7-ray ex- 
cess from point sources, but also from extended regions 
like the Orion complex which are not part of the model 
(see Fig. 3 at longitude I — 180°). At a few positions the 
emission was found to be systematically low. We decided 
to disregard these isolated regions and excluded 6 % of 
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Fig. 2. RMS deviations between observations and model cal- 
culations as described in Sect. 3. The RMS scatter derived for 
individual scans at constant longitudes is displayed as circles, 
the total RMS deviation between model and data is given by 
the crosses. 



the data, identical for all models, from the RMS determi- 
nation. 

Fig. 2 presents the derived RMS deviations for our 
models. Model 1 is not a fit but gives the residual 7-ray 
emission after subtracting the 7-ray emission from disk 
and diffuse ionized gas. For the models 2 to 5 the addi- 
tional emission originating from a baryonic halo has been 
calculated. Model 2 represents the best fit. The RMS de- 
viation between data and model (crosses) as well as the 
scatter between individual scans at constant longitudes 
(circles) is minimal. In Fig. 3 we plot model 2 in compar- 
ison with the observations. 

The models which have been represented in Fig. 2 have 
been supplemented by additional calculations with various 
core radii a and flattening parameters q according to (jfy . 
We found 0.2 £ q ^ 0.4 and a ~ 5.6 kpc to fit the obser- 
vations well, however, in no case we could recover the best 
fit results of model 2. A flattening parameter q — 0.3 cor- 
responds to the flattening of model 2, which has no free 
parameters concerning the shape of the halo. 

Halos with a flattening parameter q }t 0.4 barely fit 
the observations. In particular, spherical halo models, as 
represented by model 4 and 5, result in very poor fits. 
Flat dark matter models with q ^ 0.2 or models with a 
scale height h z ;$ 1 kpc for the CR distribution (Combes 
& Pfenniger 1996) also do not fit the EGRET data in a 
satisfactory way. 

Concerning the question, whether the 7-ray emission 
from H2 clumps suffers from obscuration as defined by /x 
in Eq. (1), we need to distinguish two cases, the transpar- 
ent and the opaque model. For model 3 we derive /x = 1- 
From Fig. 1 it is obvious that the radii of the clumps 
must be large, r c £ 20 AU. Such models have been de- 
tailed by Pfenniger & Combes (1994), by Combes & Pfen- 
niger (1996), and by Gerhard & Silk (1996). For model 2 
we derive /x = 0.09, the opaque case. The clumps have 
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Fig. 3. The EGRET diffuse 7-ray emission (E 1 > 100 MeV) 
observations - lines with dots, and best fit model 2 - solid lines. 
The longitude I varies from 0° (bottom) to 300° (top) in steps 
of AZ = 30°. 



masses of ~ 10 _3 Mq and radii r c ~ 6 AU. For masses of 
0.3 — 3 • 1O _3 M0 we derive radii between 3 and 10 AU 
respectively. The lower values are close to the radii esti- 
mated by Walker & Wardle (1998) from their explanation 
of extreme scattering events. Our upper limit corresponds 
to the estimates by Draine (1998) and by KK98. 

The major difference between model 3 and 2 is the 
assumption concerning the distribution of cosmic rays on 
large scales. Models 3 to 5 are based on a constant en- 
ergy density of 0.12 eV cm~ 3 out to distances of 100 kpc 
(De Paolis et al. 1999). In model 2 it is assumed that the 
cosmic rays are in pressure equilibrium with the observed 
gaseous halo. This results in a rather narrow distribution 
with an exponential scale height of h z ~ 4.4 kpc. Accord- 
ing to current diffusion models, the distribution of cosmic 
rays is restricted to z-scales of 2-4 kpc (Webber & Soutoul 
1998) or to 4.9±2 kpc (Ptuskin & Soutoul 1998). Similar 
parameters were used by Salati et al. (1996). Their con- 
clusion, that only about 3 % of the 7-ray flux which is 
expected for a gaseous dark matter halo can be observed 
is in good agreement with our determination of /t- Strong 



et al. (1999) favor CR re-acceleration and determine z- 
scales of 4 - 10 kpc. In this case also the inverse Compton 
emission at high latitudes would be affected. Using pa- 
rameters as proposed by Strong et al. (1999) for z-scales 
of 4 - 10 kpc leads to an increase of the RMS deviation be- 
tween model and data. Clump radii, however, are affected 
by ;$ 10% due to the steep gradient of /t (Fig. 1). 

Since there is no observational evidence for a cosmic 
ray halo at z-scales exceeding 10 kpc, we adopt our best 
fit model 2. We interpret the residual observed 7-ray emis- 
sion after subtraction of the emission from disk and dif- 
fuse ionized gas layer as due to H2 clumps with masses 
of ~ 10~ 3 M Q and characteristic radii of r c ~ 6 AU. Such 
chimps, exposed to cosmic rays, are optical thick and emit 
7-rays only close to their surfaces. The Milky Way dark 
matter halo may contain ~ 10 such H2 clumps with a to- 
tal mass of ~ 2 • 1O 11 M . The local surface column density 
of these clumps then is E ~ 140 M pc~ 2 (KK98). 
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